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Dislocation-assisted phase separation processes in binary systems subjected to irradiation effect are studied 
analytically and numerically. Irradiation is described by athermal atomic mixing in the form of ballistic flux 
with spatially correlated stochastic contribution. While studying the dynamics of domain size growth we have 
shown that the dislocation mechanism of phase decomposition delays the ordering processes. It is found that 
spatial correlations of the ballistic flux noise cause segregation of dislocation cores in the vicinity of interfaces 
effectively decreasing the interface width. A competition between regular and stochastic components of the 
ballistic flux is discussed. 
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1. Introduction 

A Study of nonequilibrium phenomena observed in materials under sustained particle or laser irra¬ 
diation attains an increasing interest in modern theoretical physics, condensed matter physics, material 
science, and metallurgy. Particle or laser irradiation causes the production of structural disorder with 
generation of a large amount of point defects. These defects can organize into defects of higher dimension 
and stimulate the occurrence of nonequilibrium phenomena. In recent decades, numerous experimental 
data have shown that alloys under sustained irradiation can be considered as nonequilibrium systems 
manifesting phase transitions, phase separation, pattern formation with rearrangement of point defects 
in bulk and on a surface (see for example, HlfTol). First observations of ordering/disordering processes in 
irradiated alloys were discussed six decades ago (see reference (ill). It was shown that nonequilibrium 
conditions in such systems are caused by interactions of high energy particles with atoms of a target (pure 
material, alloys). 

From practical viewpoint, a study of these phenomena remains an urgent problem to predict the be¬ 
havior of construction materials. A study of phase stability in various solids and metallic alloys under 
sustained irradiation received a long-standing attention due to its intrinsic interest and its relevance in 
technological problems such as: improvement of mechanical properties, radiation resistance, radiation 
damage, etc. Mechanical stability of construction materials is governed by rearrangement of the defects 
produced by irradiation and their segregation on phase interfaces and grain boundaries. Perturbation of 
the atomic configuration by irradiation causes the alternation of the phase stability 11211151 . Therefore, in 
order to predict the behavior of irradiated materials at different loading, one should know the physical 
mechanisms leading to self-organization of the defect structure that causes microstructure transforma¬ 
tions. 

It is known that phase transformations in alloys subjected to particle irradiation can be quite differ¬ 
ent from that observed in the absence of irradiation QD. Experimental observations of phase separation 
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(spinodal decomposition) in binary alloys (Ni-Cu, Ni-Cu) have shown that the electron irradiation can 
increase the solute mobility. It causes phase decomposition at temperatures at which diffusivities un¬ 
der thermal conditions are too small to provide this effect (see, for example, references (17] [18)). The 
same results were obtained for alloys Au-Ni, Cu-Ni, Fe-Mo with different compositions (T^ED. Irradia¬ 
tion damage can lead to precipitate dissolution and stagnation precipitate ordering 12211241 . It was found 
that phase separation in irradiated systems can occur at temperatures above a coherent spinodal. This 
effect can be described by point defects production with an increase of their mobility to dislocations 
in such a way that the misfit dislocations move with the composition field relieving strains. The corre¬ 
sponding model of a mobile dislocation density field coupled with the composition field was proposed in 
references 12511271 . It was found that phase separation is possible above the coherent spinodal due to the 
motion of dislocations with a decrease of misfit strains. Phase decomposition and patterning sustained 
by the dislocation field dynamics coupled with the composition field in binary systems under sustained 
irradiation were studied in reference (28). In this work, the irradiation effect was considered as an ad¬ 
ditional contribution to the free energy according to the model proposed in references 12911321 . In this 
model, the irradiation induced atomic mixing was described by a ballistic (athermal) flux responsible for 
the production of a structural disorder. In reference (28), the authors have shown that stable patterns 
characterized by time independent amplitude and wave-length emerge due to misfit dislocations. These 
linear defects are capable of reducing the coherency strain emergent at an atomic size mismatch. In the 
proposed model, the authors consider a deterministic case, where irradiation increases the free energy 
of the system. A problem related to nonequilibrium effects induced by fluctuations of the point defect 
concentration and local temperature in cascades was not solved for the system with mobile dislocations. 

In this work, we extend the above model of phase separation with a dislocation mechanism in binary 
systems subjected to irradiation taking into consideration stochastic conditions. The goal of the paper is to 
study the role of the above fluctuation effects in the prototype model of binary systems undergoing phase 
separation assisted by mobile dislocations. We take into account the stochastic component of the ballistic 
flux proposed in reference (33). Such a stochastic model was exploited to study phase decomposition 
processes 13411361 and patterning (37) in irradiated systems. The corresponding stochastic contribution 
takes care of local fluctuations in the composition field due to stochasticity of point defects concentration 
and temperature. We consider spatial correlations of these fluctuations and study the effect of spatial 
correlations onto phase decomposition processes. By taking into account the difference in time scales 
for composition and dislocation density fields we, initially, consider the simplest case of one slow mode 
(composition field). This allows us to perform the mean-field analysis for the slow mode and study the 
effect of the dislocation mechanism strength onto phase decomposition processes. The dynamics of the 
coupled, simultaneous time evolved composition and dislocation fields is studied numerically. Here, we 
discuss the statistical properties of phase separation and the domain size growth law. We show that the 
growth of domain sizes is delayed by dislocations participating in phase decomposition processes under 
sustained irradiation. This effect was predicted theoretically in unirradiated systems (see references (38] 
(39)). Considering the segregation of dislocations in the vicinity of interfaces, we discuss a competition 
between regular and stochastic components of ballistic flux fluctuations. 

The work is organized as follows. In section we present the stochastic model of a binary system 
with ordinary thermal fluctuations representing the internal noise and ballistic flux fluctuations playing 
the role of external noise. In sectionwe study a reduced model, where dislocation density is excluded 
according to an adiabatic elimination procedure. In section]^ we numerically consider the dynamics of 
the phase decomposition. We conclude in sectionj^ 


2. Model 

Considering a binary alloy A-B, one can exploit the Bragg-Williams theory, where the correspond¬ 
ing free energy density is /bw(<^a) = ZwoCaCb/ 2Tcilncf. Here, i = {A,B}, Cf = Nf/N, N = A/a + A^b 
is the total number of particles, Z is a coordination number, T is the temperature measured in ener¬ 
getic units, an ordering energy Wq - {2w{a,b} - uJ{a,a} - uJ{b,b}) is defined through pair interaction en¬ 
ergies After expanding fyw around the critical concentration c = 1/2, we arrive at Landau-like 
potential fiy/) - -Ay/^12 IA with A - r/[c(l - c)] - Zwq and B - [Zwo-l/il-c]^ l/c^]/3; the 
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quantity y/ measures the deviation from the critical concentration, i.e., y/ = ic-c). Taking into account 
the inhomogeneity of the alloy and assuming that y/ir) varies slowly on the scale of lattice parameter 
ao, i.e., y/{r-\-ao) - y/{r) + ao • Vi//(r), one can take into account the gradient energy term to the free en¬ 
ergy in the form (Vi//)^/2, where ro is the interaction radius determining the interface width between 
two phases enriched by atoms A and atoms B. Following the Krivoglaz-Clapp-Moss expression, one has 
Tq = di;(/c)/dA:^, where vik) is the Fourier transform of the atomic interaction energy. Adopting the Cahn- 
Hilliard approach, the dimensionless free energy functional assumes the Ginzburg-Landau form lioll^ 

^0 - f dV^ - ^y/^ -p -p y (Vi//)^j. The case A < 0 corresponds to temperatures above the chemical 
spinodal. 

An additional contribution to the free energy is given by a lattice mismatch in the form of elastic 
energy = {v^EI2) f dVy/^ (28). Here, E is the Young modulus, v relates to the lattice parameter change 
with respect to the composition (Vegard’s law), i.e., a = vy/) |42l|43). The elastic contribution shifts 

the corresponding coherent spinodal: A = v^E. 

Following reference (28), we take into account the dislocation-assisted mechanism for spinodal de¬ 
composition by introducing dislocation-dislocation interactions in the formed = fdV[^\b\^ ^(V^tD)^]. 

Here, the constant C relates to a core energy of dislocations. The elastic strain energy of the system is gov¬ 
erned by the Airy stress function a) satisfying the equation VtD = E{\/xby - ^ybx), where bx and by are 
the corresponding components of the continuous dislocation density field in a two dimensional problem. 
This term accounts for the nonlocal elastic interaction between dislocations. 

To describe the coupling between the composition field y/ and strain field of dislocations, we use 
the results of the work reference (28) and introduce a relevant contribution to the free energy in the 
form ^c-vf dVy/S/^(D. This two-dimensional model was previously used to study the melting (44H4^ . 
dislocation patterning (47) and phase separation in the misfitting binary thin films (48). 

By combining all the above contributions, the total free energy of the actual system reads ^tot = 
^0 + =^6 + Therefore, the dynamics of the conserved fields y/, bx and by is described by the 

following set of deterministic equations with diffusive dynamics 


oy/ 


(2.1) 

dtb:c = [MgVl + McV^y] 

1 

(2.2) 

' Sb, ’ 

dtby = [M,vl + MgVl] 

d’etat 

(2.3) 

' 6by ■ 


Here, M is the solute mobility. Mg and Me denote the mobility for glide and climb, respectivel30 

The effect of irradiation leads to an increase in the total free energy due to ballistic mixing of atoms 
in cascades. One of the models allowing one to describe these processes was proposed in reference (29). 
It is based on the introduction of the spatial coupling term relevant to ballistic exchanges under irradia¬ 
tion conditions. The related Langevin dynamics with the additive external noise mimicking a stochastic 
ballistic mixing was studied in reference (32). It should be noted that this approach does not properly 
take into account the fiuctuations of the solute by a stochastic motion of the defects in cascades. As far as 
these fiuctuations occur in a correlated medium (crystals), the corresponding spatial correlations of fiuc¬ 
tuations should be considered. The other concept of a ballistic mixing describing the above mentioned 
fiuctuations was proposed in reference 11211141 [^ . It was shown that a ballistic mixing is stochastic in 
nature since the knocked atoms move at random at the distance R. According to this approach, the bal¬ 
listic mixing can be described by introduction of the ballistic diffusion fiux with a fiuctuating ballistic 
diffusion coefficient. These fiuctuations are induced by irradiation (fiuctuations in both concentration of 
defects and local temperature in cascades). In reference (33) it was shown that such an approach leads 
to a multiplicative noise Langevin dynamics, where spatially correlated external fiuctuations promote 
the solute fiux opposite to the ordinary diffusion fiux. The phase decomposition of binary systems under 
the above assumptions was studied in references (MIISS), while patterning processes in one component 
crystalline systems under the irradiation effect were discussed in references (3^137) . 

^This set of equations belongs to the models with conserved dynamics according to the classification suggested by Galperin and 
Hohenberg in reference (49]. 
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In this paper we exploit the model of a stochastic ballistic flux according to discussions provided in 
references I12lll4l[33ti35ll37|[^ . We assume that the force mixing induced by ballistic jumps occurs with 
relocation distances b = (R) = / RiviR)dR, where R is distributed according to the known distribution 
wiR). Such ballistic jumps can be considered as a non-thermal diffusion process with a “diffusion coeffi¬ 
cient” D^. The corresponding ballistic flux is = -D^S/y/ (121. Following reference (33l, we assume that 
such a diffusion occurs in the fluctuating environment. Indeed, collision processes of an energetic parti¬ 
cle with an atom result in local fluctuations in the temperature and a number of point defects (Frenkel 
pairs). It allows one to introduce fluctuations of the ballistic flux assuming t), where 

C(r, t) is the random source. Therefore, the quantity has regular (deterministic) and stochastic contri¬ 
butions, i.e., 1b = + regular part, is characterized by the quantity Db = defined 

through a frequency of atomic jumps (pay, where cp and ay are irradiation flux and replacement cross- 
section, respectively. The corresponding stochastic part, relates to fluctuations in atomic relocation 
distances. It is characterized by a dispersion {{dR)^). Therefore, for the ballistic flux, one can write 

J^ = -[D^ + C(r,f)]Vt//. (2.4) 

Here, we assume that realizations C(r, t) are independent in time but correlated in space. Statistical prop¬ 
erties of the external noise f (r, t) are as follows: (C(r, t)} = 0, (C(r, f)C(r', t)} = a^DbCir-r')6it - t'). Here 
C(r-r') = ( v^rc)“^exp{-(r-r')^/2rc } is the spatial correlation function with the correlation radius 
re, (J^ is an external noise intensity describing a dispersion of the quantity D^. The quantity in the 
correlator (r, ^)C(r^ t)) means that external fluctuations are possible only at nonzero irradiation flux. In 
such a case, the right-hand side of equation can be written as a sum of thermally sustained diffusion 
flux J = -MV6F[y/] /6y/ and the ballistic flux J^. 

To proceed, we act onto equation l |2.2| by the operator Vy and act onto equation l |2.3| by Adding 
these two equations, we arrive at one equation for the density field cp = In our consideration, we 

take into account that the solute mobility M can depend on the field i/a as M = MoM(i//). Next, let us move 

to dimensionless quantities: y/' - \/BIAy/, a)' - yJsiEr^cD, a - y/ElAv, r' - yjAlr^r, t' - iMoA^lrQ)t, 

D\ = DbIMoA, b' = ySBr^/A^h, M'g = Mc,gr^E/MoA^, Me = Mg, e = Ad Erl, ^ = Mc,grlE! MqA^ . 
Considering a general case, let us put M[y/') = 1 -y/'^. Using the above renormalizations and dropping 
the primes, we arrive at a system of two equations 

dty/ = V 'Miy/) \dl^^niy/)Vy/-V^y/] -\-a<p-\-V •(ir,t)Vy/-\-V • y/Miy/)(ir, t), 

^ ^ (2.5) 

dt(p- -m[(p-\- aV^y/- eV^cp), 


where 




Db 

Mill/) ’ 


fW = - 


1 - 
2 







( 2 . 6 ) 


The last term in the equation for y/ represents an internal multiplicative noise. It is characterized by 
(0 = 0 and (r, t)^ (r, t)) =65 (r-r'; t- f), where 6 is the parameter measuring the internal noise intensity 
proportional to a bath temperature. We assume that no spatio-temporal correlations between fluctuation 
sources are possible. 

It should be noted that time scales of the evolution of composition and dislocation density fields de¬ 
scribed by the quantity m cx Mc,g/Mo can be different. At m = 0, we get a system with immobile disloca¬ 
tions. Limit m ^ oo corresponds to extremely mobile dislocations. It means that m e [0,oo) depends on 
the properties of the studied material and can be considered as a free parameter of the model. A detailed 
study of the systems characterized by different values of m was reported in reference (sT). Next, follow¬ 
ing reference |5T). we consider the system with mobile dislocations by taking m ^ 1. In the simplest case 
of extremely mobile dislocations (m » 1), one can adiabatically eliminate the fast field considering the 
behavior of the slow one. In our further study, we discuss statistical properties of the system according to 
subordination principle. To make a general analysis, we study the behavior of the system with the above 
two fields by taking into account the above time scales difference. 


^As far as (p is defined in terms of gradients of bx and by, we can monitor the strain energy reduction at segregation of disloca¬ 
tions at interfaces. 
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3. Subordination principle and mean-field results 


3.1. stability of the reduced system 

Let US consider the simplest case when mobile dislocations instantaneously adjust the evolving com¬ 
position field. To this end, we put m » 1. This allows us to exclude the fast variable (p by assuming 
m~^dt(p - 0. Hence, using the Fourier representation for the Fourier co mpo nents cpk and y/k, we ob¬ 
tain the relation from the second equation of the system |2.5|. In the case ek^ « 1, we 

can expand the denominator up to the first order and obtain an approximation cf)]^ ^ a k^jl + ek^)y/]c, or 
(f) ~ -aV^(l - eS/^)y/. Substituting this expression into the first equation of the system l |2.5) , we get one 
equation for the slow mode in the form 

dty/ = S/-M{xifWfi{xif)+S/-[({r,t)S/y/+ f)], (3.1) 


where the notation Vpiy/) = Vid^fiy/) - ■ 


eV^)xi/ is introduced for convenience; plays the 

role of the effective chemical potential [VjUgf = The obtained equation i3.li is the 

main equation for the reduced system analysis. According to the structure of equation l |3.1) , one should 
have in mind that for the field y/ we get conserved dynamics, i.e., /i//(r, t)dr = y/Q, where y/o stands for 
the initial concentration difference; y/o = const according to the mass conservation law. 

In statistical analysis we study only observable (averaged) quantities. By averaging equation l |3.1) 
one gets noise correlators which can be calculated using the Novikov’s theorem (52) (the correspond¬ 
ing averatging procedures are shown in appendix A). The thermal fiux (internal) noise correlator reads: 


= -(0/2)V • Calculations for the external noise correlator give: {(S/y/) = 

[C(0)V^(i//)-p C"(0)V(i//)], where we have to note that C(r-r') acquires its maximal value at r = r', 
which implies that VC(r-r')|^^^/ = 0; C(r-r')|^^^/ = C"(0) < 0 (see references (33]|34l[53][54) for de¬ 

tails). Therefore, after averaging we get 


dtixi/) = V • (MVp> - - V • (Va^M) + 0-2 [C"(0)V2{1/.) + C(0)V^i/^)]. 


(3.2) 


Let us study the stability of the homogeneous state y/-0. As far as we consider the system with con¬ 
served dynamics, the corresponding stability analysis can be done studying the dynamics of the structure 
function 5(/c, t) as the Fourier transform of the two point correlation function {y/ir, t)y/ir^, t)). Linearizing 
the system in the vicinity of the state yr = 0,m the continuous and thermodynamic limit, we arrive at the 
dynamical equation for the structure function in the form (see appendix B for details) 

dS(fc, t) ^ _2k^w{k]S{k, r) + 2gfc^ - f dqS(q, t) + f dqC(|k-q|)S(q, f), (3.3) 

at {Znr J (Znr J 

with the dispersion relation 

OJik) = Cef-^-Pefk^. (3.4) 

Here, Cgf is the effective control parameter playing the role of an effective temperature counted from the 
critical value and ^ef is the inhomogeneity parameter defined as 

£ef=a^-l + d + Db[l+a^C''{0)], Pei=l-a^e-Dba^Cm. (3.5) 


It follows that the ballistic diffusion (its regular component) increases the effective temperature of the 
system, whereas correlation effects governed by the term C"(0) < 0 decrease its value. At the same time, 
ballistic diffusion is capable of decreasing the interface width between two phases [last term in ^ef 
equation (3.5) 1. 

From equation (3.4) , one finds that the critical wave-number that bounds the unstable modes is de¬ 
fined as 


I - a^e-Dija^dO) 


(3.6) 
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Figure 1. (Color online) The structure function for the reduced system: (a) the dynamics of Sik, f) at = 
0.15, = 0.1; (b) the dependence Sik] at a fixed time interval it = 10). Other parameters are: a = 0.5, 

0 = 0.1, ^ = 0.2, rc = l. 


The most unstable mode is described by the wave-number k^n = /cq/ V2. For the actual set of the sys¬ 
tem parameters at a < 1, one gets a decreasing dependence kda). Therefore, at small a, the dislocation 
mechanism promotes a decrease in the wave-number of unstable modes. With an increase in a, spatial 
modulations of the composition field are characterized by long-wave modes. At the same time, spatial 
correlations of the external noise f increase the wave-number of unstable modes due to C"(0) < 0. 

Typical dynamics of the structure function Sik, t) are shown in figure[^(a), Ajm relates to the position 
of the peak in the dependence Sik) [see figure[^(b)]. From figure[^(a) one can see that during the system 
evolution, the peak of 5(/c, t] related to the wave-number /cm moves toward /c = 0 and its height increases. 
Therefore, the corresponding spatial instability promotes the ordering processes with the formation of 
domains of phases enriched by atoms A or B. The effect of the system parameters onto Sik) is shown 
in figure [^(b). Here, one can find that an increase in the ballistic mixing coefficient promotes the 
formation of the structural disorder characterized by realization of long-wave perturbations and small 
maximal value of the structure function. The stochastic contribution of the ballistic mixing flux acts in 
an opposite manner stimulating the ordering processes. At elevated values of external noise intensity cr^, 
the corresponding spatial structures are characterized by lower domain sizes enriched by the atoms of 
one sort. This effect is caused by spatial correlations of the external fluctuations. 

The linear stability analysis is valid only on a short time scale. At large time limit (t oo), one can 
use the mean-field approximation based on the analysis of the solution of the Fokker-Planck equation for 
the probability density of the composition field. 

3.2. Mean-field approximation 

To analytically study the statistical properties of phase separation at f ^ oo one needs to analyze a 
stationary probability density The behavior of the system can be described analytically within 

the framework of the mean-field approach derived for the systems with conserved dynamics f34l[3^[5^ 
[58). 

In the Wiess mean-field approximation, one can use the mean field (molecular field) rj = (y/) as an 
order parameter for phase transitions and phase decomposition. In such a case, one uses the transforma¬ 
tion procedure, which allows us to introduce the order parameter in d-dimensional space as follows: 

Ay/=^2di{y/)-y/). (3.7) 

The mean-field value rj is self-consistently defined according to the definition of the mean (y/) through 
the stationary distribution function P^. In the mean-field theory, the stationary distribution is a function 
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of y/ and r/. A procedure to obtain the corresponding distribution as a solution of the corresponding 
Fokker-Planck equation is shown in appendix C. 

In order to define the transition and critical points at phase separation, we use the procedure pro¬ 
posed in references (55][56]|59). According to this approach in a deterministic case with = 0, one has 
a model dty/ = V • MS/6^/6yf, where the restriction y/Q = fy dry/ir, t) is taken into account, y/o is fixed 
by the initial conditions. For such a system, the transition point is Gjiy/o): at 0 > Gjiy^o), the homoge¬ 
neous state y/Q is stable; at 0 < Gjiy/o), the system separates into two bulk phases, y/i and y/ 2 , fulfilling 
(y/) = y/Q, The transition from a homogeneous state to two-phase state is critical for y/o = 0 only, i.e. 
0t(O) = 0c is the critical point. The corresponding steady state solutions are given as solutions of the 
equation S/MS/6^I6y/ = 0. If no fiux condition is applied, then the bounded solution is 6^l6y/ = h, 
where /z is a constant effective field (in equilibrium systems h is a chemical potential). In the homoge¬ 
neous case, the value h depends on the initial conditions y/o- Above the transition point, the steady state 
is not globally homogeneous. Here, the system separates into two bulk phases with the values y/i and 
y /2 - In the case of the symmetric form of the free energy functional where two phases with y/i = -y /2 are 
realized, we get h = 0 (55). Hence, if the field h becomes trivial, then the transition point can be defined. 

By using this procedure, one finds that in the homogeneous case the mean-field is the same every¬ 
where and equals the initial value, i.e. rj = y/o- Hence, solving the self-consistency equation 


1 

r/ = y* y/Ps{y/,ri,h)dy/ 


(3.8) 


at the fixed mean-field value, we obtain the constant effective field h. Below the threshold 0 t, the system 
decomposes into two equivalent phases with (y/i) = - {y/ 2 ), and h should be the same for these two phases 
and should be zero. Hence, below the threshold only {y/} should be defined as a solution of the self- 
consistency equation with Pg iy^> 0)- 

In the actual case, we are interested in phase decomposition phenomena induced by the irradiation 
effect. Therefore, we define the transition and critical points only for the parameters relevant to irradi¬ 
ation, namely D^, by fixing y/o- The corresponding dependencies of the effective field h versus Dij 
and the external noise intensity are shown in figure From figure (a) it is seen that the field h takes 
nonzero values above the transition point Di^j. According to the definition of /z as a chemical potential, 
one can say that at fixed y/o, the quantity h is compensated and phase separation occurs inside the do¬ 
main of the parameters where h = 0. From the dependencies /z(D^) it follows that the ordered state with 
the initial concentration y/o can be found only before Dijj.lt is seen that with an increase in the strength 
of the dislocation mechanism described by a, the transition point decreases. Considering the depen¬ 
dence hia^) [see figure[^(b)], it follows that with an increase in a, the phase separation processes can be 




Figure 2. (Color online) Constant field h versus and [panels (a) and (b), respectively] at y/o = 0.05 
and different sets of the system parameters. 
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Figure 3. (Color online) The mean-field value t] versus ballistic diffusion coefficient Dij and noise intensity 
cr^ [panels (a) and panel (b), respectively]. 


realized inside the noise intensity interval 

Next, let us discuss the mean-field 77 behavior varying the system parameters. Here, we solve the 
self-consistency equation at /z = 0. According to figure [s] (a), the mean-field decreases with an increase 
in the coefficient It behaves critically in the vicinity of the value when nontrivial values 

iy^i) - -( 1 /^ 2 ) appear. Here, we arrive at the conclusion that the irradiation leads to homogenization of 
the composition field distribution. By increasing the noise intensity, one finds that phase decomposition 
is realized at lower values of compared to the case of small a^. An increase in the intensity of internal 
fluctuations 6 suppresses the phase separation at large The same effect can be found when the inten¬ 
sity of the feedback between dislocation density and composition field increases. This result follows even 
from the analysis of the deterministic system, where the elastic field changes the critical point position. 
A more interesting situation is observed by varying the external noise intensity [see figure [^(b)]. Here, 
one finds that the external noise leads to the emergence of a disordered state (77 = 0) at > a^. In other 
words, external fluctuations of large intensity lead to a statistical disorder. On the other hand, at special 
choice of the system parameters related to ballistic flux properties, a reentrant behavior of the mean-field 
rj is observed. Here, phase decomposition is realized in a window of the noise intensity e 
This phenomenon is caused by the competition between regular and stochastic (correlated) parts of the 
ballistic flux. At the most essential contribution is given by the regular component leading to 

homogenization of the alloy. Inside the interval e correlation effects dominate and lead 

to a decrease in the effective temperature of the system. At large a^, external fluctuations destroy the 
ordered state. Therefore, the correlated ballistic flux is capable of inducing phase separation processes 
of initially homogeneous alloys. According to dependencies h[a^), one can conclude that the dislocation 
mechanism sustains the above reentrance. 

The corresponding phase diagram illustrating the formation of ordered and disordered phases is 
shown in figure It is seen that a reentrant phase separation is realized in a narrow interval for 
at elevated At small Dij, one gets the standard scenario of phase decomposition, where fluctuations 
suppress the ordering processes. Comparing different curves, it follows that an increase in the intensity 
of internal fluctuations 6 shrinks the interval for the system parameters bounding the domain of the 
ordered phase (cf. dash and dash-dot-dot curves). The domain of reentrant ordering extends with an in¬ 
crease in the external noise correlation radius rc (cf. dash and dot lines). At an elevated rc, the size of the 
domain for the ordered phase grows at small An increase in a shrinks the domain of the ordered 
phase and extends the domain of the reentrant decomposition (cf. dot and dash-dot curves). 

According to the obtained results, it follows that phase separation processes can be controlled by the 
main system parameters (temperature and elastic properties of the alloy) and statistical properties of ir¬ 
radiation effect (regular and stochastic contributions in the ballistic flux). Moreover, correlated stochastic 
contribution of this flux is capable of inducing reentrant phase separation processes. 
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Figure 4. (Color online) The mean-field phase diagram at different sets of the system parameters. 

Let us study a strong coupling limit, neglecting the spatial interactions term. Assuming yz-rj, one gets 
the stationary distribution in the form Psiy/,rj) = 6[y/- rj). To obtain an equation for the effective field h, 
we integrate equation l |C.27) and find 



(3.9) 


Ath = 0, one has solutions for two bulk phases 



(3.10) 


The corresponding transition line can be obtained directly from equation ^3.10^ at 77 = y/o, where y/o is 
the initial value for the composition field. Critical values for the system parameters can be obtained from 
the condition r/ = 0. It is interesting to note that ballistic flux parameters lead to renormalization of the 
effective temperature: 0 counting off the critical one 0c = 1. Indeed, according to the definition of the 
free energy density fiy/) for an unirradiated system, the quantity 0 is reduced to a^. In an irradiated 
system, 0 is reduced to + D^-P0-2dD^c7^(Co - Q). Hence, the regular component of the ballistic flux 
increases the effective temperature in the same manner as the internal noise does. On the other hand, 
the external fluctuations reduce this temperature due to their spatial correlations. From equation l |3.10) 
it follows that an increase in D^, 6, and a causes a decrease in the order parameter. The external noise 
is capable of extending the interval for where the mean-field takes up nonzero values. 

It is known that the mean-field results are mostly qualitative. To validate the mean-field results, we 
will use a simulation procedure. In the next section we discuss the behavior of our system considering the 
dynamics of both quantities y/ and cf) and numerically illustrate a possibility of reentrant phase separation 
processes. 

4 . The effect of dislocation density field dynamics 

4 . 1 . stability analysis 

Considering the system with two fields y/ and (p, let us start with stability analysis. Averaging the 
system l |2.5) over noises, we get dynamical equations for average fields in the form 



(4.1) 
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Figure 5. Stability diagram (a) and critical wave-number dependence on and (b) at m = 1, a = 0.5, 
rc = 0.5. 


Next, let us consider the stability of the state (i// = 0, (/) = 0) using Lyapunov’s analysis for fluctuations of 
both y/ and cf). A linearization of the governing equations in the Fourier space yields 


/d^x 

dr 

m 

^ dt ' 


Ik^wik) 
\ mak^ 


-m(l -h ek^)^ 


(W 


(4.2) 


where 

wik) = a^-ld-6d-Dh[ld-a^C"iO)] -r [l-D^a^C(O)] k^. (4.3) 

The corresponding Lyapunov exponent takes the form 


2 A= [/c^ir'(/c)-m(l-h^/c^)] ± [k‘^w{k) ek'^)Y dm [l -r k'^wik) k^]. (4.4) 

According to the analysis of the Lyapunov exponent j4.4| , we can And critical values for and a^, 
bounding the domain of unstable modes. The corresponding stability diagram is shown in flgure[^(a). It 
is seen that spatial instability characterized by A(/c) > 0 is possible only inside the window for the ballistic 
diffusion coefficient. At the same time, a growth in shrinks the interval for where A(/c) > 0. From 
the stability diagram, one finds that at small the external noise can sustain a spatial instability even at 
large intensities of fluctuations. An increase in the correlation radius rc of these fluctuations enlarges the 
instability domain. Therefore, strongly correlated external fluctuations are capable of inducing spatial 
instability at short time scales. As figure [^(b) shows, the critical wave-number bounding wave-number 
of unstable modes decreases with and cr^. Therefore, at a large ballistic mixing intensity and the 
intensity of external fluctuations, long-wave spatial instabilities should emerge over the whole system. 

Let us consider the behavior of the structure function 5(/c, t). To obtain a dynamical equation for 
5(/c, t) in the vicinity of the point {y/ -^,(p - 0), we exploit the approach previously described by con¬ 
sidering the system of two equations l |2.5) . Moving to the discrete representation and multiplying every 
equation from the system by y/, we arrive at the system of two equations 


dSikyt) T r T 2A:^ 

—= - 2k^w(k)S(k, t) + aG(k, t) + 2ek^ -^ 

d? {27iY 

dGik, t) 


j dqSiq, 


t) 


2k^Dha^ 

{27iY 


j dqC(|k-q|)S(q,?), 


dt 


= -2m [(l -r ek^) G(/c, t) - ak^Sik, t )], 


(4.5) 


where G(/c, t) = {y/]^{t)(p-]^{t)) - {y/-]^{t)(p\^{t)), and w{k) is given by equation l |4.3| . 


23003-10 












































Phase separation processes in binary systems subjected to irradiation 




Figure 6. (Color online) The structure function behavior for the original system: a) the dynamics of Sik, t) 
at Dij = 0.15, (7^ = 0.1; b) the dependence S{k) at f = 10 and a different set of the system parameters. 
Other parameters are: 6 = 0.1, a = 0.5, e = 0.2, Tq = 1.0, m = 1.0. 


The dynamics of Sik, t) is shown in figure (a). Comparing graphs for 5(/c, t) related to actual and 
reduced models [cf. figures (a), (a)], one finds that the peak of the structure function is larger in the 

actual two-component dynamical model. From the obtained dependencies for the structure function Sik) 
shown in figure [^(b) it is seen that an increase in the regular component of the ballistic flux essentially 
decreases the structure function; it shifts the peak position toward small wave-numbers. Considering 
the effect of the dislocation density mechanism strength, one finds that with an increase in a, the wave- 
number of unstable modes decreases [see the insertion in figure [^(b)]. At the same time, the height of 
the peak of the structure function decreases at a short time scale. Therefore, the dislocation mechanism 
is capable of delaying the ordering processes. Herein below, we will show that this effect can be observed 
by the dynamics of the average domain size. 

4.2. Numerical results 

To qualitatively describe the system behavior, we numerically solve the system l [2.5) . In simulation 
procedure, the Heun method was used. The system was studied on the lattice with square symmetry of 
the linear size L = 128^ with periodic boundary conditions and the mesh size £ = 0.5; At = 0.005 is the 
time step. We take (i//(0,0)) = ((/)(0,0)) = 0, ((5i//(0,0))^) = ((50(0,0))^) = 0.01 as initial conditions. The 
obtained results are statistically independent of different realizations of noise terms f, f. 

Typical evolution of both the composition field y/ and dislocation density cf) is shown in figurej^ Here, 
the regions of high values of both fields y/ and 0 are represented by white, whereas the black areas 
relate to lower values of the corresponding field. The coupling between dislocation density and compo¬ 
sition field is well observed. A coordinated motion of dislocations and phase boundaries was previously 
observed at atomic scale using phase field models (see references l28l[60l ) 

It is seen that the dislocation field takes up large values in the vicinity of the phase boundaries; inside 
the decomposed phases, the dislocation density is around zero. The corresponding oscillations of y/ near 
the interfaces indicate that the strain energy is reduced due to the atomic size mismatch (28). Therefore, 
misfit dislocations segregate on the boundaries. In figurej^ we plot the oscillating structure of dislocation 
density field cf) corresponding to the distribution of the concentration field. It is seen that in the vicinity 
of the boundary, (p changes the sign; inside the phases, <p^0. 

To make a quantitative analysis of phase separation, let us study the dependencies of dispersions of 
the fields y/ and cp defined as Jy/ = (iSy/)^} and = (idcp)^}, where dy/ -y/- (y/), d(p - (p- {(p). At {y/) - 
{(p) = 0, these quantities are reduced to the second statistical moments playing the role of effective order 
parameters at phase decomposition processes (due to conservation laws for y/ and (p). The quantity Jy/it) 
is proportional to the area below the structure function Sj^it), i.e., Jy/it) = Therefore, the growth 
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t=10 t=100 t=1000 t= 10000 



Figure 7. Evolution of both the composition field y/ (top panels) and dislocation field (bottom panels) at 
Dij = 0.15, 0-2 = 0.5, 6 = 0.1, a = 0.5, ^ = 0.2, = 1.0, m = 1. 


in (and the related growth of 7^) corresponds to the phase decomposition process. Both moments 
and J(p grow toward nonzero stationary values. The related stationary values ^ oo) can 

be used to define two phases separated during the the long time evolution of the system. If dislocations 
are excluded from the description, then one can use only to monitor the formation of two phases. In 
our case, the dynamics of phase decomposition can be described by an additional order parameter 
manifesting segregation of dislocations in the vicinity of the interface. From naive consideration, one can 
expect that at large values of both and 7^, one gets two well decomposed phases with an increased 
dislocation density at the interface. At small and 7^, one gets a mixed state where no large deviations 
in the composition field are observed. On the other hand, it means that dislocations are distributed inside 
the phases. 

In figurej^ we plot the dependencies of the order parameters at different values of the external noise 
intensity at other fixed system parameters. It is seen that an increase in the noise intensity results in 
small values of both and 7^. At = 0 (see solid lines in figure]^, the order parameter grows toward 
stationary value in a monotonous manner. This means an increase in the area under the corresponding 
structure function 5(/c, t) and the formation of well decomposed phases enriched by atoms of different 
sorts (see right-hand panel illustrating the distribution of the composition field y/). The order parameter 
J(p initially grows meaning the formation of two separated phases with an increasing dislocation density 
in the vicinity of the interface. At the next stage, a decaying dependence of 7^ is observed. This means 



Figure 8. (Color online) Snapshots of both concentration and dislocation fields at t= 10000 shown in 3D- 
view, illustrating the oscillating behavior of dislocation density field near the boundaries of two phases 
y/ = +l and y/ = -1. Other parameters are: Djj = 0.15, = 0.1 ,6 = 0.01, a = 0.5, e = 0.2, Vc = 1.0, m = 10. 
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Figure 9. (Color online) The dynamics of both order parameters Jy/ and at different values of the ex¬ 
ternal noise intensity and the corresponding snapshots of both concentration and dislocation density 
fields att = 4000. Other parameters are: Djj = 0.15, 6 = 0.01, a = 0.5, e = 0.2, Tq = 2.0, m = 10. 


an agglomeration of the domains belonging to one phase resulting in annihilations of dislocations with 
opposite signs. At the late stage, the dislocation density goes to its stationary value together with Jy;. At 
a small noise intensity (see dashed curves in figurej^, one observes the same dynamics of both order 
parameters, where Jy^ and take up low values. Here, the external noise sustains the formation of an 
ordered state characterized by separated phases. This effect is caused by correlation properties of the 
external noise. The deterministic part of the ballistic flux acts in the manner opposite to the stochastic 
contribution. By increasing the noise intensity (in the domain of a disordered phase according to the 
mean-field analysis), one gets a transition toward disordered state. Here the order parameter Jyr attains 
a very small stationary value (see dotted curves in figure [^. The formation of a disordered state is well 
accompanied by time independence of the quantity 7^. It fastly attains a small stationary value and fluc¬ 
tuates around it. Here, there are no phases enriched by atoms A or B (see the right-hand panel in figurej^, 
dislocations are distributed over the whole system. Therefore, the effect of fluctuations characterized by 
large values of becomes larger than the correlation effects, leading to the formation of a totally disor¬ 
dered state. 

Let us consider in detail the effect of the dislocation density field onto the dynamics of the phase de¬ 
composition. An evolution of both order parameters Jy; and at different values of the coupling constant 
a is shown in figure [To| (a). Here, one can see that both Jy; and increase with a (the order parameter 
J(p has the corresponding peak at transition to the coarsening regime). Let us consider the behavior of 
the stationary order parameter - Jy^it ^oo) [see the insertion in the upper panel in figure [l^ (a)]. It 
rapidly increases at small a and slowly grows with a. From the obtained results it follows that a strong 
coupling between the composition field and the dislocation density field urges the formation of the or¬ 
dered state due to redistribution of dislocations over the whole system, as well as their motion to the in¬ 
terfaces. Therefore, phase separation is well sustained by a dislocation field. It is interesting to compare 
the dynamics of the average domain size at different a. According to discussions provided in references 
(38][39l it is known that dislocation mechanism is capable of changing the dynamical exponent z, describ¬ 
ing the domain size growth law (R) oc To analyze the dependence we calculate the averaged 

value {k{t)) oc l/{R{t)} according to the standard definition (/c(f)) = //c5(/c, t)dk/ JS{k)dk. In the stan¬ 
dard theory of phase decomposition, the corresponding Lifshitz-Slyozov approach gives z = 1/3 (61). The 
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(b) 


Figure 10. (Color online) (a) Dynamics of both order parameters Jy/ and at different a and Dij = 0.12, 
(7^ = 0.01. (b) Evolution of the average domain size at different a and = 0.1, = 0.5. Other parame¬ 

ters are: 0 = 0.01, e = 0.2, Tq = 2.0, m = 10. 


same value for z is observed when phase separation is sustained by vacancy mechanism. If dislocation 
mechanism of phase decomposition plays the major role, then the dynamical exponent takes up lower 
values z - 1/6 (39). In our case, we can control the strength of the dislocation mechanism varying param¬ 
eter a. According to the results in figure [To|(b), one obtains z - 0.33 at a ^ 0, as Lifshitz-Slyozov theory 
predicts. This result was obtained for the system subjected to a stochastic ballistic flux with another form 
of the function M(i//) (see reference (34l)- It was shown that an increase in the external noise intensity 
results in disordering processes. Comparing the curves related to a = 0.1 and a = 0.5, it follows that the 
dislocation mechanism delays the dynamics of the domain sizes growth. In our case, at a = 0.5, we get a 
lower value for the dynamical exponent. 

Finally, let us consider the dynamics of the average dislocation density {(pn) Id the vicinity of the in¬ 
terfaces and the coherence (interface) width (L^), where this density decreases toward zero value inside 
the separated phases. We calculate the quantity {(pn) as the mean height of (p{r) profile averaged over 
the whole system. The coherence width (L^) is calculated as the width of the interface on the half-height 
of (p[r) profile averaged over the system. In figure pT] (a), solid and dashed lines denote one-dimensional 
profiles for the composition and dislocation density fields, respectively. The dynamics of both {(pn) and 
{Lcp) is shown in figure [ll|(b). Comparing the data related to different sets of and one finds that 
the dislocation density attains a stationary value during the decomposition process. Considering the dy¬ 
namics of {(pn) at different external noise intensities, it follows that the growth in cr^ increases the values 
of the dislocation density. Therefore, the external noise is capable of inducing a phase decomposition ac¬ 
companied by segregation of dislocations at interfaces. On the other hand, with an increase in D^, the 
quantity {(pn) takes up lower values. Here, dislocations are distributed over the whole system due to 
homogenization of the system produced by a regular component of the ballistic flux. The competition be¬ 
tween regular and stochastic components of the ballistic flux can be observed by studying the dynamics 
of the averaged interface width. From the bottom panel in figure (b), one finds that during the sys¬ 
tem evolution, the quantity (L^) attains a maximum. This maximum relates to the stage of the domains 
growth. A coalescence regime (large domains absorb small ones) is accompanied here by an increase 
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Figure 11. (Color online) (a) Illustration of typical one-dimensional profiles for solute concentration y/ 
and the field (p taken at the time t = 400 = 0.12, = 0.1). (b) Evolution of both the average maxi¬ 

mum dislocation density field {(pn) a coherence length (L^p) at different values of and Other 
parameters are: a = 0.5, e = 0.2, rc = l, m = l,6 = 0.1. 


in the interface width. A decrease in (L^) corresponds to a coarsening regime. At large time intervals, 
{L(p} attains a stationary value. Comparing the corresponding stationary values at different it follows 

that a regular component of the ballistic flux produces a disordering accompanied by extension of the 
interface width (cf. the curves marked by circles and triangles). Comparing the curves with different 
one finds that the noise is capable of playing a constructive role leading to a decrease in the interface 
width. Localization of dislocations at interfaces can be stimulated by a correlation effect of the external 
fluctuations [see the curve with triangles in the upper panel in figure [TT](b)]. It should be noted that the 
above effect is possible only inside the domain of the ordered phase formation. Therefore, the mentioned 
above constructive role of external fluctuations is caused by their spatial correlations. 


5. Conclusions 

We have studied phase separation processes driven by the dislocation evolution mechanism in a bi¬ 
nary system subjected to sustained irradiation. We describe the irradiation effect by introducing a ballis¬ 
tic flux having stochastic properties. We assumed that these fluctuations are spatially correlated. 

By taking into account different time scales of evolution of both composition and dislocation density 
fields, we have initially considered a reduced model using the adiabatic elimination procedure. In this 
case, the dynamics of the composition field playing the role of a slow mode is studied. By making use of 
the linear stability analysis, we have shown the constructive role of ballistic flux fluctuations. These fluc¬ 
tuations induce spatial instability at a short time scale. At a large time scale, we have used the mean-field 
approach allowing us to describe the properties of phase separation. Corresponding mean-field phase di¬ 
agrams illustrating the possibility of phase decomposition are calculated. It was found that ballistic flux 
components reduced to the intensity of atomic mixing (ballistic diffusion coefficient) and the intensity 
of the corresponding fluctuations are capable of controling reentrant phase separation processes. It was 
shown that a reentrant character of phase separation is governed by spatial correlations of the external 
noise and mobile dislocations. 

Considering the dynamics of both studied fields by means of computer simulations, we have found 
that the formation of domains enriched by atoms of different sort is accompanied by an increase of the 
dislocation density field in the vicinity of the interface. Inside such domains, dislocation density takes up 
zero values and all dislocations characterized by different signs in two phases segregate on interfaces. 
The average length of the interface decreases at the formation of the ordered state. In the disordered 
state, all dislocations are distributed over the whole system. It was shown that spatially correlated ex¬ 
ternal fluctuations act in the manner opposite to the regular component of the ballistic flux and induce 
the ordered state formation. Studying the effect of the dislocation density field onto phase decomposition 
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we have considered the dynamics of the average domain size at different strengths of the dislocation 
mechanisms. It was shown that the universal dynamics of the average domain size delays due to a re¬ 
distribution of dislocations. We have found that at small contribution of the dislocation density field, the 
corresponding universal dynamics is described by the standard Lifshitz-Slyozov law with dynamical ex¬ 
ponent z = 0.33. With an increase in the dislocation mechanism strength, this exponent takes up lower 
values and the domains of new phases are characterized by smaller linear sizes. 
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Appendix A 


Let us represent the system on a regular d-dimension lattice. Within the framework of the standard 
formalism of a discrete representation, the system can be divided onto cells of the linear size L = 
£N, where £ is the a mesh size. Then, the partial differential equation US is reduced to a set of usual 
differential equations written for every cell i on a grid in the form 

^ (A.1) 

where the index i labels cells, i = \ .the discrete left-hand and right-hand operators are: 




(A.2) 


Discrete correlators of stochastic sources are of the form: 


where Cf-j is the discrete representation of the spatial correlation function C(r) which in the limit of 
zero correlation length becomes 6ij/£^. For the two-dimensional problem considered below, the quantity 
C\i-j\ can be computed as a discrete version of the Fourier transform of C(r-r') written in the form 
(53) C(/c) = exp{-(rc/2)[sin^(/Cx/2) -r sin^(/Cy/2)]}. Noise correlators can be calculated according to the 
recipes shown in references (34l|53]|54||58l. Next, we consider the case d = 2. 

To calculate the internal noise correlator ( we use the Novikov theorem which can be written 

in the form 



Using a relation 


5^/Wi 




dy/Wi 


dipiit) 


t=t' 


and a formal solution of the Langevin equation (A.l), one can write 


dy/iit) 






Substituting it into equation l|A.4), we get 


Mi^t) } = 


£^ \ dy/j 


dA/M; Q /dM 

'Mj ) = — ' ■' 


2^2 >J \ dy/j / • 


(A.4) 


(A.5) 
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Using the relation between left-hand and right-hand discrete gradient operators and moving to the con¬ 
tinuum limit, we get 

= (A.6) 

To calculate the external noise correlator {y/Oy we again use the Novikov theorem written as follows: 


I 

f ((jWM 


Sy/jit) 

s(j{n 


dt' ^Dba^Y.Cj-1 


/ Sy/jit) 


(A.7) 


According to the formal solution of the Langevin equation, the corresponding derivative with respect to 
f takes the form {6y/j/6C]c\t=f) - ^jk Next, substituting it into equation l |A.7| and using a discrete 
representation of the Laplacian, we find the sum over the index k allowing us to write 

[Ciiy^j+i) + -2Co{i/^j)]. (A.8) 

Acting by Laplacian operator onto this construction, we finally obtain in continuum limit 

My/O = DfeO-2 + C"(0)A{i/^)]. (A.9) 


Appendix B 


To obtain a dynamical equation for the structure function, we need to construct a dynamical equa¬ 
tion for the two-point correlation function {y/iy/j}. In our computation procedure, we use the procedure 
described in reference (62). Multiplying the linearized Langevin equation l |A.l| written for y/j onto y/f 
and doing the same procedure for the linearized Langevin equation written for y/f, we add these two 
equations. This procedure allows us to obtain a dynamical equation for {y/iy/j} written in the form 

^ik[[Db -1 + a^) iyfjyrk) - (1 - ^klW0j) + Wjy^k^k)] 
{yfi'^)k^k)-'2-{'iyi'^)0k^k) 

^jk [[Db-i+iyfiyfk) - (i - ^kiWiyfi) + (y'iy'k(k)] 

{yrjV\k^k) -2{yfjV]^yrk(k). (B.IO) 


diy/jy/j) 

dt 

-P 

-P 

-P 


where the sum runs over the repeating indexes. Using the Novikov theorem with recipe shown in Ap¬ 
pendix A, one can calculate the following correlators: 


6 6 

(yfi^)k^k) (yfi^)0k^k) = [^]k(y'ky'j)+^ik{y'ky'j)), 

(y^jfk(k) = Dba^ [AkiCk-iiy/iy'j) + AjiCk-iiyriy/k)) ■ 


Introducing the structure function 5(k, t) = {y/kit)y/-kit)) in the discrete space 

Sy{t) = {N£y^{y/y{t)y/-y{t)) 

with Fourier components 

yrAt) ^ y/mit) = 


one can define the derivative 


dSv ik rr-r 

—- = {£!N)^^ —-— 
dt dt 


In the following computations, one needs to exploit definitions: 

£ \2 

- e^^^^^J-^^^Ais{y/sy^j) = DySyit), 


(B.12) 
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(B.13) 

2(^1 e‘'''''’-^-'-'>C|,-_y|{i/A,-VAy) = (JV/)-2^C^S^_v(r), (B.14) 

where Dy = £~^['L\i\=i cos(kvri) - l) can be understood as the Fourier transform of the discrete Lapla- 
cian. Using the above definitions and equation jB.lO) with correlators we arrive at the equation 

2Dy — 1 + 0 + DIj + DIjCT^ /^QjjiCjji — Dy Sy{t) 

2eDy + 2Dym)-^J^S^it)-2DhaHN£)-^DyJ^Cy.^S^it). (B.15) 

In continuous and thermodynamical limit (N ^ 00 , £ ^ 0), we obtain a dynamical equation for the struc¬ 
ture function in the form 

-2k^ {a^ - I6-\-Dha-^(J^C"m-\-k^ [I-a^e-Dha^dO)]} 5(/c, t) 

2ek^--^^ f dqS(q, t) -h f dqC(|k-q|)5(q, r). (B.16) 

(ZTTj^ J (ZTTj^ J 


d5(/c, t) 
dt 


for the quantity Sy{t): 

dSyit) 

dt 


Appendix C 


To obtain the probability density distribution ^{[y/], t) in d-dimensional space, let us start with defi¬ 
nitions: 

j \ _ 

m{nf\,t) = (Upi(t)\ = {p(t)), Pi(t)=5(yi(t)-yi\^, (C.17) 

where .. .ic and (...) are the averages over initial conditions and fiuctuations, respectively. To obtain the 
corresponding Fokker-Planck equation, we use a standard technique and exploit the stochastic Liouville 
equation for the distribution pit) in the form 

d . 

(3rp = --— iy/ip). (C.18) 

oy/i 

Inserting the expression for the time derivative from equation \A.l) and averaging over the noise, we get 

^ ■ (C.19) 

Correlators in the second term can be calculated by means of the Novikov theorem, that at ^ = 1 
gives (52) 


(v]j y^Kj{t)p)=e j dt'sj,5{t-1') 


<A 


L 

ijy/j(j(t)p}^Di,a^ j dt'C\j-ic\Sit- 


s(kit') /■ 


f) 


Introducing notations gij - A = {f.O, for the last multiplier, one has 

SgijPit) _ d dyriit) I 

SAkit') ~ y^'hyk8Ak(t')\t=/' 


Using a formal solution of the Langevin equation, the response functions take up the form 


5nti(t) 


t-t' 


5nti(t) 

S(kit') 


t=t' 


(C.20) 


(C.21) 


(C.22) 
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After some algebra, we obtain the Fokker-Planck equation for the total distribution ^ in the discrete 
space 


dyr, 


dif/j d-iffi ^ dy/i 


where the relations between left-hand and right-hand gradient operators are used. 

To proceed, let us obtain an evolution equation for the single-site probability distribution 


(C.23) 


After integration one has 


where 


Piit) = [ \ 




dPjit) d 
dt dy/, 






df + Dij 1 , „ _ , 

-- y/j -p —- 1- 

A/T. 2d\ M/ ' 


diifj Mj 

d 


2 


6 \/Mj ^ \/Mj -P DijG ^mn^lj-nlV^n- 


^ dy/n 


In the stationary case with no flux, one arrives at the equation 

Aij{.y^j)Psiy/i) == 0, 


(C.24) 


(C.25) 


(C.26) 


where Pg is a stationary distribution function. By taking i - j, dropping the subscripts, we arrive at the 
mean-fleld stationary equation EH) 


■hP, 


M 


M 


2 


^ ^ a^-pD^ , „ _ , 


- -dy/M -\-2dDija^y/ 


M 


CiT]- -Co—1// 

dy/ dy/ 


Ps. 


(C.27) 


Here we have used the mean-fleld approximation, allowing us to write 


^0^7 “I ^ V^nn{i) 2dy/i 

\nn{i) 


■2di{xl/}-xl/), r]=(-if/}, 


(C.28) 


where nn{i) denotes the nearest neighbors of a given site. The mean-fleld value rj and the integration 
constant h should be deflned self-consistently. Equation jC.27) has the solution in the form 


Ps [y/j ri,h)-N exp 


,0(i/A',r/, h) 


r 


(C.29) 


where 


Q 

0(i/A,r/, h) - -Mdy;f - (a^ -pD^)i//-p [M-a^e] {rj-y/)- -dy^M -ZdO^a^Coy/ h, 
Eiy/; rj) = 6M -P 2dDija^y/iCoy/ - Cirj) 

the normalization constant is a function of the mean-fleld rj and the constant h: 


(C.30) 


N = N{rj, h)^ 


^ , \ 
,0(i//',r/, h) 




(C.31) 
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Aoc/iiA>KeHH5i npoL^eciB (t>a30Boro posiuapyBaHHB 3a HaBBHOCTi 
Ai/ic/ioKai^iM B 6iHapHi/ix CMCTeiviax, niAA^HMX onpoMiHCHHio 

fl.O. XapHCHKO, O.M. moKOTOBa, A.L BaiuTOBa. 1.0. /li/iceHKO 

iHCTHTyT npMKTiaAHoT c|)i3MKi/i HAH VKpaiHi/i, By;i. nerponaByiiBCbKa 58, 40000 CyMi/i, VKpaiHa 

ripoBeAeHo Aoc;iiA>KeHH?i npoi^eciB (t)a30Boro po3LuapyBaHH?i 3a Ai/ic;ioKai4iCiHi/iM MexaniaMOM b 6iHapHi/ix ci/i- 
CTeMax, nifl^aHi/ix am onpoMiHeHH^i. OnpoMiHeHH^i onucycTbca arepiviiMHi/iM nepeMiiuyBaHH^iM aroMiB, aa pa- 
xyHOK yBeAeHHfi 6a/iicTi/iMHoro noToxy, u^o Mae npocTopoBO-CKope/ibOBaHy CTOxacTi/iMHy CK/iaAOBy. flpi/i bh- 
BMeHHi Ai/iHaMiKM pocryAOMeHiB noKaaaHO, u^o Awc/iOKai^iMHUM MexaniaM ynoBmbHKDe npoi^ecynop^iAKyBaHH?!. 
BcTaHOB/ieHO, u^o npocTopoBi Kope/i^M^ii LuyMy 6a/iicTi/iMHoro noToxy CTi/iMy/iK)K)Tb cerperai^iK) 9\fl,ep Ai^c/ioxa- 
u}\A B oxoyii Mix(t)a3Hnx rpaHni^b, e(t)exTi/iBHo 3MeHLuyK)MM Lui/ipi/iHy Mix(t)a3Horo Luapy. Poar/iJiMyro xoHxypeH- 
i^iK) Mix peryyifipHOK) ra cToxacTUMHOK) xoMnoHeHiaMM 6a;iicTi/iMHoro noToxy. 

K/iiOHOBi c/ioea: (pasoBe posLuapyBaHHfi, onpoMiHeHHfi, inyM 
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